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Abstract. We update gamma-ray burst (GRB) luminosity relations among certain 
spectral and light-curve features with 139 GRBs. The distance modulus of 82 GRBs at 
z > 1.4 can be calibrated with the sample at z < 1.4 by using the cubic spline interpolation 
method from the Union2.1 Type la supernovae (SNe la) set. We investigate the joint 
constraints on the Cardassian expansion model and dark energy with 580 Union2.1 SNe la 
sample (z < 1.4) and 82 calibrated GRBs data (1.4 < z < 8.2). In ACDM, we find that 
adding 82 high-z GRBs to 580 SNe la significantly improves the constrain on Vt m — 
plane. In the Cardassian expansion model, the best fit is Vt m = 0.24lon5 and n = 0.16^o; 3 2 
(lcr), which is consistent with the ACDM cosmology (n = 0) in the 1 a confidence region. 
We also discuss two dark energy models in which the equation of state w(z) is parametrized 
as w(z) = w 0 and w(z) = w 0 + W\z/{ 1 + z), respectively. Based on our analysis, we see 
that our Universe at higher redshift up to z = 8.2 is consistent with the concordance model 
within lcr confidence level. 

1 Introduction 

In recent years, the combined observations of nearby and distant Type la supernovae (SNe 
la) have provided strong evidence for the current accelerated expansion of the universe 
EM- The cause of the acceleration remains unknown. Many authors suggest that the 
composition of the Universe may consist of an extra component called dark energy, which 
may explain the acceleration of the Universe at the current epoch. For example, the dark 
energy model with a constant equation of state Pjp = w = — 1 is one of the several possible 
explanations for the acceleration. While other models suggest that dark energy changes with 
time, and there are many ways to characterize the time variation of dark-energy. Here, we 
adopt a simple model in which the dark-energy equation of state can be parameterized by 
P/p = w(z) = wo + W\z/{1 + z) = tco + wi(l —a) EE- where wo is constant, w\ represents 
the time dependence of dark energy, and a = 1/(1 + z) is the scale factor. In addition, 
models where general relativity is modified can also drive universe acceleration, such as the 
Cardassian expansion model is a possible alternative for explaining the acceleration of the 
universe that invokes no vacuum energy [6]. 

SNe la have been considered a perfect standard candle to measure the geometry and 
dynamics of the Universe. Unfortunately, the farthest SNe la detected so far is only at 
z = 1.914 [Tj. It is difficult to observe SNe at z > 2, even with excellent space-based 
platforms such as SNAP jS]. And this is quite limiting because much of the most interesting 
evolution of the Universe occurred well before this epoch. Gamma-ray bursts (GRBs) are the 
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most luminous transient events at cosmological distances. Owing to their high luminosities, 
GRBs can be detected out to very high redshifts [2j. In fact, the farthest burst detected so 
far is GRB 090423, which is at z = 8.2 [TO]. 0 Moreover, in contrast to SNe la, gamma-ray 
photons from GRBs are almost immune to dust extinction, so the observed gamma-ray flux 
is a direct measurement of the prompt emission energy. Hence, GRBs are potentially a more 
promising standard candles than SNe la at higher redshifts. The possible use of GRBs as 
cosmological probes started to become reality after some empirical luminosity relations were 
discovered. These GRB luminosity relations have been proposed as distance indicators, such 
as the correlations ri ag — L [12], V — L [13], E p — E iso [H], E p — L [T5l [16], E p — E P [17] . 
Trt — L [IS] , and so on. Here the time lag (ri ag ) is the time shift between the hard and soft 
light curves; the luminosity (L) is the isotropic peak luminosity of a GRB; the variability 
(V) of a burst denotes whether its light curve is spiky or smooth, and V can be obtained by 
calculating the normalized variance of an observed light curve around a smoothed version 
of that light curve [13]; ( E p ) is the burst frame peak energy in the GRB spectrum; (Ei SO ) 
is the isotropic equivalent gamma-ray energy; (E y ) is the collimation-corrected gamma-ray 
energy; and the minimum rise time (trt) in the gamma-ray light curve is the shortest time 
over which the light curve rises by half of the peak flux of the pulse. However, Ref. [T9] 
found that the updated V — L correlation was quite scattered. Its intrinsic scatter has been 
larger than the one that could be expected of a linear relation. 

Generally speaking, with these luminosity indicators, one can make use of them as stan¬ 
dard candles for cosmological research. For example, Ref. [2I[J constructed the first GRB 
Hubble diagram based on nine GRBs using two GRB luminosity indicators. With the 
E p — E 1 relation, Ref. [21] placed tight constraints on cosmological parameters and dark 
energy. Ref. [22] used a model-independent multivariable GRB luminosity indicator to 
constrain cosmological parameters and the transition redshift. Ref. [18] made use of five lu¬ 
minosity indicators calibrated with 69 events by assuming two adopted cosmological models 
to construct the GRB Hubble diagram. Ref. [23] suggested that the time variation of the 
dark energy is small or zero up to z ~ 6 using the E p — L relation. Ref. [24] extended the 
Hubble diagram up to z = 5.6 using 63 gamma-ray bursts (GRBs) via E p — L relation and 
found that these GRB data were consistent with the concordance model within 2a level. 
In a word, a lot of other works in this so-called GRB cosmology field have been published 
(please see [19] and [25] for reviews). However, there is a so-called circularity problem in 
the calibration of these luminosity relations. Because of the current poor information on 
low-z GRBs, these luminosity relations necessarily depend on the assumed cosmology. Some 
authors attempted to circumvent the circularity problem by using a less model-dependent 
approach, such as the scatter method [26] [27], the luminosity distance method [25] , the 
Bayesian method [23 . 50] . and the method by fitting relation parameters of GRBs and cos¬ 
mological parameters simultaneously [3T] |32]. However, these statistical approaches still can 
not avoid the circularity problem completely, because a particular cosmology model is re¬ 
quired in doing the joint fitting. This means that the parameters of the calibrated relations 
are still coupled to the cosmological parameters derived from a given cosmological model. 

To solve the circularity problem completely, one should calibrate the GRB relations in a 
cosmology independent way. Recently, a new method to calibrate GRBs in a cosmological 
model-independent way has been presented [33] |34j, 35]. This method is very similar to the 
calibration for SNe la by measuring Cepheid variables in the same galaxy, and it is free from 

1 A photometric redshift of 9.4 for GRB 090429B was reported by ;TTi. 
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the circularity problem. Cepheid variables have been regarded as the first order standard 
candles for calibrating SNe la which are the secondary standard candles. Similarly, if we 
regard SNe la as the first order standard candles, we can also calibrate GRBs relations with 
a large number of SNe la since objects at the same redshift should have the same luminosity 
distance in any cosmology. This method is one of the interpolation procedures which obtain 
the distance moduli of GRBs in the redshift range of SNe la by interpolating from SNe la 
data in the Hubble diagram. Then, if we assume that the GRB luminosity relations do 
not evolve with reshift, we can extend the calibrated luminosity relations to high -2 and 
derive the distance moduli of high- z GRBs. From these obtained distance modulus, we can 
constrain the cosmological parameters. 

In this paper, we will try to determine the cosmological parameters and dark energy 
using both the updated 139 GRBs and 580 SNe la. In Section 2, we will describe the 
data we will use and our method of calibration. To avoid any assumption on cosmological 
models, we will use the distance moduli of 580 SNe la from the Union2.1 sample to calibrate 
five GRB luminosity relations in the redshift range of SNe la sample (z < 1.4). Then, the 
distance moduli of 82 high- z GRBs (z > 1.4) can be obtained from the five calibrated GRB 
luminosity relations. The joint constraints on the Cardassian expansion model and dark 
energy with 580 SNe and 82 calibrated GRBs data whose z > 1.4 will be presented in 
Section 3. Finally, we will summarize our findings and present a brief discussion. 

2 Calibrating the updated luminosity relations of GRBs 

2.1 Observational data and methodology 

As mentioned above, we calibrate the updated luminosity relations of GRBs using low -z 
events whose distance moduli can be obtained by those of Type la supernovae. Actually, 
we use the cosmology-independent calibration method developed by Refs. [33) [MJ [35] • This 
method is one of the interpolation procedures which use the abundant SNe la sample to 
interpolate the distance moduli of GRBs in the redshift range of SNe la sample (z < 1.4). 
More recently, the Supernova Cosmology Project collaboration released their latest SNe la 
dataset known as the Union2.1 sample, which contains 580 SNe detections [36]. Obviously, 
there are rich SNe la data points, and we can make a better interpolation by using this 
dataset. 

Our updated GRB sample includes 139 GRBs will redshift measurements, there are 57 
GRBs at z < 1.4 and 82 GRBs at z > 1.4. This sample is shown in Table [Q which includes 
the following information for each GRB: (1) its name; (2) the redshift; (3) the bolometric 
peak flux P b oic>; (4) the bolometric fluence Sboioj (5) the beaming factor /beam! (6) the time 
bg Ti ag ; (7) the spectral peak energy E p ; and (8) the minimum rise time trt- All of these 
data were obtained from previously published studies. Before GRB 060607, we take all the 
data directly from Ref. [18]. We adopt the data between GRB 060707 and GRB 080721 
from Ref. [T9J . For those GRBs detected after July 7th, 2008, we adopt the data directly 
from Ref. [37]. Applying the interpolation method, we can derive the distance moduli of 57 
low-z GRBs and calibrate five GRB luminosity relations with this low- z sample, i.e., the 
Pag — L relation, the E p — L relation, the E p — P 7 relation, the trt — L relation, and the 
E p — E iso relation. The isotropic peak luminosity of a burst is calculated by 

L = 4nDlP holo , (1) 
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the isotropic equivalent gamma-ray energy is given by 

E iso = 4nDlS holo (l + z) -1 , 

and the collimation-corrected energy is 

-®7 -^iso/beam 47T J D^ bolo / beam (l + z)- 1 . (3) 

Here, is the luminosity distance of the burst, P bo io and Shoio are the bolometric peak flux 
and fluence of gamma-rays, respectively, while / beam = (1 — cosset) is the beaming factor, 
djet is the jet half-opening angle. We assume each GRB has bipolar jets, and P 7 is the true 
energy of the bipolar jets. 

For convenience, the luminosity relations involved in this paper can be generally written 
in the power-law forms 

log y = a + b log x, (4) 

where a and b are the intercept and slope of the relation, respectively; y is the luminosity 
(L in units of erg s^ 1 ) or energy (Pise or P 7 in units of erg); x is the GRB parameters 
measured in the rest frame, e.g., Ti ag (l + ^)^ 1 /(0.1 s), P p (l + z)/(300 keV), P p (l + ^)/(300 
keV), trt(1 + z) -1 /(0.1 s), P p (l + 2 ;)/ (300 keV), for the 5 two-variable relations above. 


Table 1: Luminosities and luminosity indicators. 


GRB 

z 

^bolo 

(erg cm - ^ s — ■*■) 

*5bolo 
(erg cm - 

/ beam 

T lag 

(s) 

Ep 

(keV) 

r RT 

(s) 

970228 

0.70 

7.30E-06 

± 

4.30E-07 









115± i 

0.26 

± 

0.04 

970508 

0.84 

3.30E-06 

± 

3.30E-07 

8.09E-06 

± 

8.10E-07 

0.0795 

± 

0.0204 

0.50 

± 0.30 

389± i 

0.71 

± 

0.06 

970828 

0.96 

1.00E-05 

± 

1.10E-06 

1.23E-04 

± 

1.20E-05 

0.0053 

± 

0.0014 



298±|R 

0.26 

± 

0.07 

971214 

3.42 

7.50E-07 

± 

2.40E-08 







0.03 

± 0.03 

190±|g 

0.05 

± 

0.02 

980613 

1.10 

3.00E-07 

± 

8.30E-08 









92±g 




980703 

0.97 

1.20E-06 

± 

3.60E-08 

2.83E-05 

± 

2.90E-06 

0.0184 

± 

0.0027 

0.40 

± 0.10 

254± Ia 

3.60 

± 

0.50 

990123 

1.61 

1.30E-05 

± 

5.00E-07 

3.11E-04 

± 

3.10E-05 

0.0024 

± 

0.0007 

0.16 

± 0.03 

604± l 




990506 

1.31 

1.10E-05 

± 

1.50E-07 







0.04 

± 0.02 

283±® 

0.17 

± 

0.03 

990510 

1.62 

3.30E-06 

± 

1.20E-07 

2.85E-05 

± 

2.90E-06 

0.0021 

± 

0.0003 

0.03 

± 0.01 

126±« 

0.14 

± 

0.02 

990705 

0.84 

6.60E-06 

± 

2.60E-07 

1.34E-04 

± 

1.50E-05 

0.0035 

± 

0.0010 



189±j| 

0.05 

± 

0.02 

990712 

0.43 

3.50E-06 

± 

2.90E-07 

1.19E-05 

± 

6.20E-07 

0.0136 

± 

0.0034 



65±m 




991208 

0.71 

2.10E-05 

± 

2.10E-06 









190± 1fi 

0.32 

± 

0.04 

991216 

1.02 

4.10E-05 

± 

3.80E-07 

2.48E-04 

± 

2.50E-05 

0.0030 

± 

0.0009 

0.03 

± 0.01 

318±™ 

0.08 

± 

0.02 

000131 

4.50 

7.30E-07 

± 

8.30E-08 









163± 

0.12 

± 

0.06 

000210 

0.85 

2.00E-05 

± 

2.10E-06 









4 °8±« 

0.38 

± 

0.06 

000911 

1.06 

1.90E-05 

± 

1.90E-06 









9 86±J™ 

0.05 

± 

0.02 

000926 

2.07 

2.90E-06 

± 

2.90E-07 









lOOitf 

0.05 

± 

0.03 

010222 

1.48 

2.30E-05 

± 

7.20E-07 

2.45E-04 

± 

9.10E-06 

0.0014 

± 

0.0001 



309±^ 

0.12 

± 

0.03 

010921 

0.45 

1.80E-06 

± 

1.60E-07 







0.90 

± 0.30 

89±ip 

3.90 

± 

0.50 

011211 

2.14 

9.20E-08 

± 

9.30E-09 

9.20E-06 

± 

9.50E-07 

0.0044 

± 

0.0011 



59± l 




020124 

3.20 

6.10E-07 

± 

1.00E-07 

1.14E-05 

± 

1.10E-06 

0.0039 

± 

0.0010 

0.08 

± 0.05 


0.25 

± 

0.05 

020405 

0.70 

7.40E-06 

± 

3.10E-07 

1.10E-04 

± 

2.10E-06 

0.0060 

± 

0.0020 



364±™ 

0.45 

± 

0.08 

020813 

1.25 

3.80E-06 

± 

2.60E-07 

1.59E-04 

± 

2.90E-06 

0.0012 

± 

0.0003 

0.16 

± 0.04 

142± I! 

0.82 

± 

0.10 

020903 

0.25 

3.40E-08 

± 

8.80E-09 









2- 6 ±f! 




021004 

2.32 

2.30E-07 

± 

5.50E-08 

3.61E-06 

± 

8.60E-07 

0.0109 

± 

0.0027 

0.60 

± 0.40 

80±§j 

0.35 

± 

0.15 

021211 

1.01 

2.30E-06 

± 

1.70E-07 







0.32 

± 0.04 

46± I 

0.33 

± 

0.05 

030115 

2.50 

3.20E-07 

± 

5.10E-08 







0.40 

± 0.20 

83±* 

1.47 

± 

0.50 

030226 

1.98 

2.60E-07 

± 

4.70E-08 

8.33E-06 

± 

9.80E-07 

0.0034 

± 

0.0008 

0.30 

± 0.30 

97± S 

0.70 

± 

0.20 

030323 

3.37 

1.20E-07 

± 

6.00E-08 









44 ±i 

1.00 

± 

0.50 

030328 

1.52 

1.60E-06 

± 

1.10E-07 

6.14E-05 

± 

2.40E-06 

0.0020 

± 

0.0005 

0.20 

± 0.20 

126±R 




030329 

0.17 

2.00E-05 

± 

1.00E-06 

2.31E-04 

± 

2.00E-06 

0.0049 

± 

0.0009 

0.14 

± 0.04 

68±|;f 

0.66 

± 

0.08 

030429 

2.66 

2.00E-07 

± 

5.40E-08 

1.13E-06 

± 

1.90E-07 

0.0060 

± 

0.0029 



35±] 2 

0.90 

± 

0.20 

030528 

0.78 

1.60E-07 

± 

3.20E-08 







12.50 ± 0.50 

32±f 

0.77 

± 

0.20 

040924 

0.86 

2.60E-06 

± 

2.80E-07 







0.30 ± 0.04 

67±| 

0.17 

± 

0.02 

041006 

0.71 

2.50E-06 

± 

1.40E-07 

1.75E-05 

± 

1.80E-06 

0.0012 

± 

0.0003 



63±P 

0.65 

± 

0.16 

050126 

1.29 

1.10E-07 

± 

1.30E-08 







2.10 ± 0.30 

47+ 8 

3.90 

± 

0.80 

050318 

1.44 

5.20E-07 

± 

6.30E-08 

3.46E-06 

± 

3.50E-07 

0.0020 

± 

0.0006 



*?±l l 

0.38 

± 

0.05 

050319 

3.24 

2.30E-07 

± 

3.60E-08 










0.19 

± 

0.04 

050401 

2.90 

2.10E-06 

± 

2.20E-07 







0.10 ± 0.06 

118±P 

0.03 

± 

0.01 

050406 

2.44 

4.20E-08 

± 

1.10E-08 







0.64 ± 0.40 


0.50 

± 

0.30 

050408 

1.24 

1.10E-06 

± 

2.10E-07 







0.25 

± 0.10 


0.25 

± 

0.08 

050416 

0.65 

5.30E-07 

± 

8.50E-08 









!5±|:| 

0.51 

± 

0.30 

050502 

3.79 

4.30E-07 

± 

1.20E-07 







0.20 ± 0.20 

93±|| 

0.40 

± 

0.20 
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050505 

4.27 

3.20E-07 

± 

5.40E-08 

6.20E-06 ± 8.50E-07 

0.0014 

± 0.0007 




70_i_ 24 

81±l 4 

344±“ 

0.40 ± 0.15 

050525 

0.61 

5.20E-06 

± 

7.20E-08 

2.59E-05 

± 

1.30E-06 

0.0025 ± 

0.0010 

0.11 

± 

0.02 

0.32 ± 0.03 

050603 

2.82 

9.70E-06 

± 

6.00E-07 







0.03 

± 

0.03 

0.17 ± 0.02 

050802 

1.71 

5.00E-07 

± 

7.30E-08 










0.80 ± 0.20 

050820 

2.61 

3.30E-07 

± 

5.20E-08 







0.70 

± 

0.30 

246±« 

2.00 ± 0.50 

050824 

0.83 

9.30E-08 

± 

3.80E-08 










11.00 ± 2.00 

050904 

6.29 

2.50E-07 

± 

3.50E-08 

2.00E-05 

± 

2.00E-06 

0.0097 

± 

0.0024 




436±2po 

41±J 

498*1 

0.60 

± 

0.20 

050908 

3.35 

9.80E-08 

± 

1.50E-08 










1.50 

± 

0.30 

050922 

2.20 

2.00E-06 

± 

7.30E-08 







0.06 

± 

0.02 

0.13 

± 

0.02 

051022 

0.80 

1.10E-05 

± 

8.70E-07 

3.40E-04 

± 

1.20E-05 

0.0029 

± 

0.0001 




sio±|o 

161±?io 

0.19 

± 

0.04 

051109 

2.35 

7.80E-07 

± 

9.70E-08 










1.30 

± 

0.40 

051111 

1.55 

3.90E-07 

± 

5.80E-08 







1.02 

± 

0.10 

3.20 

± 

1.00 

060108 

2.03 

1.10E-07 

± 

1.10E-07 










65: ^6p0 

62 ±fa 

139±So 

237±“ 

0.40 

± 

0.20 

060115 

3.53 

1.30E-07 

± 

1.60E-08 










0.40 

± 

0.20 

060116 

6.60 

2.00E-07 

± 

1.10E-07 










1.30 

± 

0.50 

060124 

2.30 

1.10E-06 

± 

1.20E-07 

3.37E-05 

± 

3.40E-06 

0.0021 

± 

0.0002 

0.08 

± 

0.04 

0.30 

± 

0.10 

060206 

4.05 

4.40E-07 

± 

1.90E-08 







0.10 

± 

0.10 

75± li 

““it 00 

230±op 

156 ±fo 

95±|° 

25±| 

1.25 

± 

0.25 

060210 

3.91 

5.50E-07 

± 

2.20E-08 

1.94E-05 

± 

1.20E-06 

0.0005 

± 

0.0001 

0.13 

± 

0.08 

0.50 

± 

0.20 

060223 

4.41 

2.10E-07 

± 

3.70E-08 







0.38 

± 

0.10 

0.50 

± 

0.10 

060418 

1.49 

1.50E-06 

± 

5.90E-08 







0.26 

± 

0.06 

0.32 

± 

0.08 

060502 

1.51 

3.70E-07 

± 

1.60E-07 







3.50 

± 

0.50 

3.10 

± 

0.30 

060510 

4.90 

1.00E-07 

± 

1.70E-08 













060526 

3.21 

2.40E-07 

± 

3.30E-08 

1.17E-06 

± 

1.70E-07 

0.0034 

± 

0.0014 

0.13 

± 

0.03 

0.20 

± 

0.05 

060604 

2.68 

9.00E-08 

± 

1.60E-08 







5.00 

± 

1.00 

40± ifb 

169±3° 0 

120±” 0 

63±® 

103±|° 

0.60 

± 

0.20 

060605 

3.80 

1.20E-07 

± 

5.50E-08 







5.00 

± 

3.00 

2.00 

± 

0.50 

060607 

3.08 

2.70E-07 

± 

8.10E-08 







2.00 

± 

0.50 

2.00 

± 

0.20 

060707 

3.43 

1.53E-07 

± 

2.12E-08 

3.41E-06 

± 

1.96E-07 










060714 

2.71 

2.30E-07 

± 

1.42E-08 

6.88E-06 

± 

2.47E-07 










060729 

0.54 

1.93E-07 

± 

1.30E-08 

6.43E-06 

± 

3.16E-07 







61± 




060814 

0.84 

1.83E-06 

± 

4.44E-08 

4.94E-05 

± 

4.91E-07 




0.29 

± 

0.03 

257±|| 

80± 7JO 

451±?f 2 

2°±“ 

1.65 

± 

0.24 

060904B 

0.70 

4.37E-07 

± 

2.28E-08 

4.05E-06 

± 

2.17E-07 




0.36 

± 

0.09 

1.00 

± 

0.16 

060908 

2.43 

6.69E-07 

± 

3.36E-08 

7.68E-06 

± 

1.85E-07 




0.26 

± 

0.06 

0.52 

± 

0.09 

060926 

3.21 

1.56E-07 

± 

1.22E-08 

5.47E-07 

± 

3.80E-08 




1.03 

± 

0.11 




060927 

5.60 

4.02E-07 

± 

1.54E-08 

2.37E-06 

± 

8.67E-08 




0.12 

± 

0.04 

72+' 
399±}i 

90+^ 

51 7±|| 
606±|| 

0.46 

± 

0.12 

061007 

1.26 

7.20E-06 

± 

1.11E-07 

2.24E-04 

± 

1.72E-06 




0.11 

± 

0.01 

0.38 

± 

0.02 

061110A 

0.76 

9.79E-08 

± 

1.35E-08 

2.71E-06 

± 

1.18E-07 










061110B 

3.44 

1.79E-07 

± 

2.66E-08 

6.12E-06 

± 

3.38E-07 




0.24 

± 

0.36 

0.79 

± 

0.64 

061121 

1.31 

8.04E-06 

± 

1.07E-07 

6.53E-05 

± 

5.76E-07 




0.03 

± 

0.01 

0.98 

± 

0.19 

061222B 

3.36 

2.29E-07 

± 

3.15E-08 

5.01E-06 

± 

2.49E-07 







49±| 




070110 

2.35 

1.12E-07 

± 

1.36E-08 

4.04E-06 

± 

1.64E-07 







HO+30 




070208 

1.17 

1.39E-07 

± 

2.06E-08 

1.06E-06 

± 

1.46E-07 










070318 

0.84 

4.10E-07 

± 

2.12E-08 

7.34E-06 

± 

2.01E-07 







454±}9 

83 ±U 

31± C 

233±; 

0.72 

± 

0.24 

070411 

2.95 

1.50E-07 

± 

1.31E-08 

6.29E-06 

± 

2.19E-07 










070506 

2.31 

1.67E-07 

± 

1.38E-08 

5.16E-07 

± 

3.43E-08 




2.52 

± 

0.04 

0.12 

± 

0.06 

070508 

0.82 

7.67E-06 

± 

1.18E-07 

7.26E-05 

± 

6.15E-07 




0.04 

± 

0.01 

0.20 

± 

0.01 

070521 

0.55 

2.09E-06 

± 

5.26E-08 

2.97E-05 

± 

4.00E-07 




0.04 

± 

0.01 

222±}^ 

180 ±J 

11204- 

70+P 

44±| 

0.58 

± 

0.06 

070529 

2.50 

3.32E-07 

± 

5.08E-08 

7.44E-06 

± 

4.31E-07 










070611 

2.04 

1.45E-07 

± 

2.25E-08 

9.52E-07 

± 

8.44E-08 










070612A 

0.62 

2.77E-07 

± 

4.24E-08 

2.72E-05 

± 

9.37E-07 







2.49 

± 

1.48 

070714B 

0.92 

3.24E-06 

± 

1.46E-07 

8.91E-06 

± 

6.77E-07 




0.03 

± 

0.01 

0.45 

± 

0.04 

070802 

2.45 

6.38E-08 

± 

9.69E-09 

6.50E-07 

± 

7.05E-08 










070810A 

2.17 

2.77E-07 

± 

1.77E-08 

1.59E-06 

± 

8.43E-08 




1.09 

± 

0.23 

0.73 

± 

0.22 

071003 

1.10 

4.71E-06 

± 

1.82E-07 

6.73E-05 

± 

1.48E-06 




0.38 

± 

0.05 

799±« 

0.88 

± 

0.07 

071010A 

0.98 

1.17E-07 

± 

2.67E-08 

4.97E-07 

± 

6.05E-08 







27-1-10 

52±| 

24±| 




071010B 

0.95 

9.20E-07 

± 

2.18E-08 

8.37E-06 

± 

1.16E-07 




0.84 

± 

0.04 

1.21 

± 

0.03 

071031 

2.69 

7.08E-08 

± 

8.61E-09 

2.19E-06 

± 

1.92E-07 










071117 

1.33 

2.71E-06 

± 

5.83E-08 

7.97E-06 

± 

2.02E-07 




0.60 

± 

0.01 

'' 3=1 -3Q 

73±}j 

28±| 

0.20 

± 

0.02 

071122 

1.14 

6.76E-08 

± 

2.06E-08 

1.41E-06 

± 

1.63E-07 










080210 

2.64 

2.57E-07 

± 

1.95E-08 

4.17E-06 

± 

1.41E-07 




0.53 

± 

0.17 

0.57 

± 

0.44 

080310 

2.43 

1.83E-07 

± 

1.72E-08 

5.49E-06 

± 

2.90E-07 







0.41 

± 

0.55 

080319B 

0.94 

1.55E-05 

± 

1.91E-07 

5.25E-04 

± 

3.94E-06 




0.02 

± 

0.01 

651+1 

307+22 

0.14 

± 

0.01 

080319C 

1.95 

2.22E-06 

± 

7.79E-08 

1.77E-05 

± 

2.99E-07 







0.21 

± 

0.12 

080330 

1.51 

1.33E-07 

± 

1.80E-08 

8.77E-07 

± 

1.26E-07 







2o± l 

259±of 

!7° ±i 

7 o_)_ It) 
80±« 
66±| 
42*5 

85± S 




080411 

1.03 

1.04E-05 

± 

1.31E-07 

8.75E-05 

± 

2.01E-07 




0.21 

± 

0.01 

0.65 

± 

0.01 

080413A 

2.43 

1.22E-06 

± 

2.65E-08 

9.86E-06 

± 

1.71E-07 




0.13 

± 

0.03 

0.23 

± 

0.03 

080413B 

1.10 

3.17E-06 

± 

8.25E-08 

8.00E-06 

± 

1.52E-07 




0.23 

± 

0.01 

0.50 

± 

0.03 

080430 

0.77 

4.60E-07 

± 

2.15E-08 

3.01E-06 

± 

1.53E-07 




0.68 

± 

0.08 

0.76 

± 

0.12 

080516 

3.20 

2.77E-07 

± 

2.80E-08 

5.88E-07 

± 

5.50E-08 




0.15 

± 

0.01 




080520 

1.55 

8.23E-08 

± 

1.00E-08 

1.59E-07 

± 

3.00E-08 










080603B 

2.69 

7.57E-07 

± 

2.63E-08 

7.02E-06 

± 

1.78E-07 




0.08 

± 

0.01 

0.22 

± 

0.03 

080605 

1.64 

5.99E-06 

± 

1.10E-07 

4.72E-05 

± 

4.32E-07 




0.11 

± 

0.01 

246± 14 

394 ±J 

73±|f 

485 ±|? 

31 3 c: - 1 - 73 . 6 
OlO.Oltyp e 

93.1±ff;J 

0.22 

± 

0.01 

080607 

3.04 

8.35E-06 

± 

2.42E-07 

1.00E-04 ± 

0.00E+00 




0.04 

± 

0.01 

0.18 

± 

0.06 

080707 

1.23 

1.68E-07 

± 

1.02E-08 

1.26E-06 

± 8.87E-08 










080721 

2.60 

9.57E-06 

± 

5.01E-07 

5.99E-05 

± 3.04E-06 




0.13 

± 

0.05 

0.09 

± 

0.04 

080810 

080913 

3.35 

6.70 

9.76E-07 

2.31E-07 

± 

± 

8.40E-08 

4.00E-08 

1.80E-05 

1.15E-06 

± 

± 

1.30E-06 

1.20E-07 










080916A 

0.69 

1.06E-06 

± 

1.60E-07 

2.19E-05 

± 7.30E-06 







109+J 

248±H 

134±| 

451+22 

416±| 




081121 

2.51 

1.87E-06 

± 

4.80E-07 

1.73E-05 

± 3.20E-06 










081222 

2.77 

2.62E-06 

± 

3.80E-07 

1.70E-05 

± 

1.40E-06 










090102 

1.55 

3.90E-06 

± 

5.60E-07 

3.66E-05 

± 3.40E-06 










090323 

3.57 

3.80E-06 

± 

6.90E-07 

1.48E-04 ± 2.00E-05 
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090328 

0.74 

7. 

090423 

8.20 

3. 

090424 

0.54 

1. 

090516 

4.11 

5. 

090618 

0.54 

8. 

090715B 

3.00 

8. 

090812 

2.45 

2. 

090902B 

1.82 

2. 

090926 

2.11 

1. 

090926B 

1.24 

5. 

091018 

0.97 

1. 

091020 

1.71 

1. 

091029 

2.75 

2. 

091127 

0.49 

4. 

091208B 

1.06 

3. 


28E-06 ± 5.20E-07 
24E-07 ± 6.80E-08 
80E-05 ± 1.20E-06 
93E-07 ± 3.40E-08 
58E-06 ± 9.60E-07 
96E-07 ± 2.49E-07 
17E-06 ± 2.40E-07 
72E-05 ± 3.00E-07 
82E-05 ± 3.00E-07 
55E-07 ± 7.90E-08 
82E-06 ± 8.40E-07 
63E-06 ± 2.00E-07 
82E-07 ± 1.60E-08 
03E-06 ± 1.70E-07 
47E-06 ± 6.00E-07 


1.37E-04 ± 8.00E-06 
1.17E-06 ± 3.20E-07 
5.85E-05 ± 2.10E-06 
1.09E-05 ± 1.60E-06 
3.39E-04 ± 2.50E-05 
1.09E-05 ± 1.60E-06 
3.15E-05 ± 4.10E-06 
3.78E-04 ± 3.00E-06 
1.51E-04 ± 7.70E-06 
1.66E-05 ± 5.00E-07 
4.33E-06 ± 7.60E-07 
1.68E-05 ± 3.70E-06 
5.84E-06 ± 4.70E-07 
2.65E-05 ± 5.00E-07 
7.78E-06 ± 8.80E-07 


653 ±- 

82 ±lf 

177±j[ 

185.6±g ; | 

155.5±J?f 

134 ±|?' 

572 ±^ 

775±7} 

314± t 

78.3±f 
19.2+j* 
47.9±p \ 
6i.4±i?:| 
36±l 
124+-4 
i/4± 2Q ,l 


2.2 Calibration 


First of all, we obtain the distance moduli of 57 low-z (z < 1.4) GRBs by using cubic spline 
interpolation from the 580 Union2.1 SNe la compiled in Ref. [36]. The interpolated distance 
moduli fi of these 57 GRBs and their corresponding errors are shown in Fig. QJa). The 
SNe la data are also plotted in Fig. [T{a) for comparison. When the cubic spline interpolation 
is used, the error of the interpolated distance modulus /j for the GRB at redshift z can be 
calculated by 


a 




z i+ i 


Zi+1 


e l,i + 


Zi+1 


"H,i +1 


(5) 


where and e^+i are errors of the SNe at nearby redshifts Zi and z i+ 1 , respectively. With 
D l in units of Mpc, the predicted distance modulus is defined as 


1 + = 5 log(T> L ) + 25. (6) 

Having estimated the distance moduli /i of 57 low- z GRBs in a model independent way, we 
can convert /i into luminosity distance D l by using Eq. ([6]). From Eqs. ()T|)-()3j) with the 
corresponding Pboio, S'boio, and /beam, we can calculate L, Ei SO , and E 1 . In Fig. HXb)-[T|(f ), 
with the interpolation results, we show the hve luminosity indicators for these 57 GRBs 
at z < 1.4. For each relation, we perform a linear least-squares fit, taking into account 
both the X axis error and the Y axis error. We also measure the scatter of each relation 
with the distance of the data points from the best-fit line, as done by Ref. [38] . The best¬ 
fitting results of the intercept a and the slope b with their la uncertainties and the linear 
correlation coefficients for each relation are summarized in Table [2] The best-fitting results 
derived by using the interpolation method are carried out with these 57 GRBs at z < 1.4. 
In other word, the results are derived by using data from 27, 55, 12, 40, and 42 GRBs for 
the Ti ag — L, E p — L , E p — £7., r RT — L, and E p — E iso relations, respectively. 

Ref. [TO] found no statistically significant evidence for the redshift evolution of the lumi¬ 
nosity relations. If the GRB luminosity relations indeed do not evolve with redshift, we can 
extend the calibrated luminosity relations to high- z {z > 1.4) and derive the luminosity (L) 
or energy (Pi SO or P 7 ) of each burst at liigh-z by utilizing the calibrated relations. Therefore, 
the luminosity distance Z+ can be derived from Eqs. (JTj)-(]3j). The uncertainty of the value 
of the luminosity or energy deduced from each relation is 
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log [Ep(l+z) / 300 keV] 





Figure 1: (a): The Hubble diagram of 580 SNe la (red triangles) and 57 GRBs at z < 1.4 (black 
dots) whose distance moduli are derived by using cubic spline interpolation, (b)-(f): The T\ ag — L, 
E p — L, E p — E. y, trt — L, and E p — E[ so correlations. The five correlations are calibrated with the 
sample at z < 1.4 using cubic spline interpolation. The solid lines show the best-fitting results, 
while the dashed lines represent their 2 a dispersion around the best fits. 
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Table 2: Best-fitting results 


Relation 


a 



b 


N 

r 

Rag 

- L 

52.32 

± 

0.03 

-1.10 

± 

0.05 

27 

-0.75 


-L 

52.08 

± 

0.02 

1.60 

± 

0.04 

55 

0.79 

E P ~ 

- Ej 

50.57 

± 

0.01 

1.76 

± 

0.03 

12 

0.95 

VRT 

-L 

52.85 

± 

0.03 

-1.51 

± 

0.06 

40 

-0.66 

Ep- 

E- 

-^lSO 

52.81 

± 

0.02 

1.53 

± 

0.04 

42 

0.73 


where cr a , <7b, and a x are la uncertainties of the intercept a, the slope b, and the GRB 
measurable parameters x, and cr int is the systematic error in the fitting that accounts for the 
extra scatter of the luminosity relations. Then, we obtain the distance moduli /j for these 
82 GRBs at z > 1.4 using Eq. (JEJ). The propagated uncertainties will depend on whether 
Pboio or S'boio are given by 


a 


A 4 


2 aiogi ) 


2 


+ 


f 5 ^olo 

V21n 10 Pboio 


21 1/2 


( 8 ) 


or 


and 


^ = 


:<J\o g E k 


2 


+ 


( 5 ff gbclo 

V 2 In 10 S bolo 


21 1/2 


J 




/5 \ 2 

^2°" lo g E 7 J 


( 5 ^SbouA 2 . 

V 2 In 10 S holo J 


( 5 C/beam \ 

V 2 In 10 /beam/ 


21 1/2 


(9) 


( 10 ) 


Here we ignore the uncertainty of z in our calculations. 

After obtaining the distance modulus of each GRB using one of these relations, we use 
the same method as Ref. [TBj to calculate the real distance modulus, which is the weighted 
average of all available distance modulus. The real distance modulus for each burst is 


/^fit 


Ya~ 2 


(ii) 


with its corresponding uncertainty cr /tfit = c’’/. 2 ) -1 ^ 2 , where the summation runs from 

1 to 5 over the relations with available data, //,; and cr Mi are the best estimated distance 
modulus and its corresponding uncertainty from the zth relation. 

Fig. [2] shows the Hubble diagram from the Union2.1 SNe la sample and 139 GRBs. The 
combined Hubble diagram is consistent with the concordance cosmology. The 57 GRBs at 
z < 1.4 are obtained using interpolation method directly from SNe data. The 82 GRBs at 
z > 1.4 are obtained by utilizing the five relations calibrated with the sample at z < 1.4 
using the cubic spline interpolation method. 


3 Constraints from supernovae and GRBs 

The latest Type la SNe dataset known as the Union2.1 sample was recently released by 
the Supernova Cosmology Project collaboration, which contains 580 SNe detections (see 
















52 



z 


Figure 2: Hubble diagram of 580 SNIa (red triangles) and 139 GRBs (dots) obtained using the 
interpolation method. The 57 GRBs at z < 1.4 are obtained by interpolating from SNe la data 
(black dots); and the 82 GRBs at z > 1.4 (blue dots) are obtained with the five correlations 
calibrated with the sample at z < 1.4 using the cubic spline interpolation method. The vertical 
dotted line represents z = 1.4. The solid curve represents the best-fit cosmology for a flat ACDM 
universe: il rn = 0.29, Ha = 0.71. 


El). With luminosity distance Z)l(£, z) in units of Mpc (where £ stands for all the cosmo¬ 
logical parameters that define the fitted model), the theoretical distance modulus /i t h can 
be calculated by using Eq. (]6j) . The likelihood functions can be determined from the y 2 
statistic, 


XSN 


e 


N 


E 


[/^th(^i) Robs(^i)] 

+ ^ext + Sample ’ 


( 12 ) 


where Ci c is the propagated error from the covariance matrix of the lightcurve fit, and /i 0 b s is 
the observational distance modulus. The uncertainties due to host galaxy peculiar velocities, 
Galactic extinction corrections, and gravitational lensing are included in (j ext , and cr sample is 
a floating dispersion term containing sample-dependent systematic errors. The confidence 
regions can be found through marginalizing the likelihood functions over Hubble constant 
H 0 (i.e., integrating the probability density p oc exp(—y 2 /2) for all values of H 0 ). 

Gamma-ray bursts (GRBs) are the most luminous transient events in the cosmos. Owing 
to their high luminosity, GRBs can be detected out to the edge of the visible Universe, 
constituting a powerful tool for constructing a Hubble diagram at high- z. We use the above 
calibration results obtained by using the interpolation methods directly from SNe la data. 
The x 2 value for the 82 GRBs at z > 1.4 is given by 


2 _ [/Ah {Zi) h-fit,/] 2 

Xgrb - ~2 

4=1 vz 

where and o u . are the fitted distance modulus and its error for each burst. We also 
marginalize the nuisance parameter H 0 . 

Motivated by these significant updates in the observations of SNe la and GRBs, it is 
natural to consider the joint constraints on cosmological parameters and dark energy with 


(13) 
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the latest observational data. We combine SNe la and GRBs by multiplying the likelihood 
functions. The total y 2 value is 

Xtotal = XsNe + XgRB' (14) 

The best-fitting values of cosmological model are obtained by minimizing y 2 otal . 


3.1 ACDM model 

In a Friedmann-Robertson-Walker (FRW) cosmology with mass density U m and vacuum 
energy density Da, the luminosity distance is given as 


D l (z) = 


(1 + z) 


Hn 


sinn 


D; 


, 1/2 


d^ 


X 


D; 


0 O?77,(1 + z)3 + Da + D^.(l + Z)‘ 


(15) 


where c is the speed of light, Hq is the Hubble constant at the present time, D^ = 1 — D m —Ha 
represents the spatial curvature of the Universe, and sinn is sinh when D& > 0 and sin when 
D& < 0. For a flat Universe with D& = 0, Eq. (fT5j) simplifies to the form (1 + z)c/Hq times 
the integral. In this ACDM model, the transition redshift satisfies 


z T 


2Da\ 1/3 
- Kn s 


- 1. 


(16) 


We use the data sets discussed above to constrain cosmological parameters. In the left 
panel of Fig. |3l we show the confidence regions for (D m , Da) from 82 GRBs (dark cyan 
dash-dotted lines), 580 SNe la (blue dotted lines), and 82 GRBs + 580 SNe la (red solid 
lines), respectively. We can see that adding 82 high redshift GRBs (z > 1.4) to 580 SNe 
la (z < 1.4) significantly improves the constraint on D m — Da plane. The 1 a confidence 
region from all the data sets is (D m ,DA) = (0.27to;oe, 0.62to;Jg) with y 2 lin = 727.01 for 659 
degrees of freedom. Under the assumption of a flat universe (solid line), the contours yield 
(D m ,DA) = (0.29jZo!o 4 , 0.71±q;q 4). The transition redshift at which the Universe switched 
from deceleration to acceleration phase is Zt = 0.64^q) ) | at the la confidence level (the 
right panel of Fig. |3]). 


3.2 Cardassian Expansion Model 

Ref. [6] proposed the Cardassian expansion model as a possible alternative for explaining 
the acceleration of the Universe that invokes no vacuum energy. This model allows an 
acceleration in a flat, matter-dominated cosmology. If we consider a spatially flat FRW 
Universe, the Friedmann equation is modified as 

V 2 = ^(p + Cp”). (17) 

This modification may arise as a consequence of embedding our observable Universe as a 
(3+1) dimensional brane in extra dimensions or the self-interaction of dark matter. The 
luminosity distance in this model is given by 

D l {z) = cHq\ 1 + z) r dz[{ 1 + z) 3 n m + (1 - D m )(l + z) 3n }~ 1/2 . (18) 

Jo 
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Figure 3: Left: The contour confidence levels of (f2 m , Da) in the ACDM model from the data for 
82 GRBs (z > 1.4) (dark cyan dash-dotted lines), 580 SNe la (blue dashed lines), and 82 GRBs 
+ 580 SNe la (red solid lines), respectively. The contours correspond to 1, 2, and 3 a confidence 
regions. Right: The probability versus the transition redshift derived from the GRB and SNe la 
sample. 


Fig. [4] shows constraints on Vt m and n from la to 3cr confidence regions by fitting obser¬ 
vational data. The dark cyan dash-dotted lines and blue dashed lines represent the results 
from 82 GRBs and 580 SNe la, respectively. The red solid contours show the constraints 
from the combination of these data. The best values are VL m = 0.24ldii and n = 0.161^52 
at the la confidence level with Xmin = 727.31 for 659 degrees of freedom. This result is 
consistent with the ACDM cosmology (n = 0) in the la confidence region. 

3.3 w(z) = ivq Model: Constant Equation of State 

For the dark energy model with a constant equation of state (w(z) = Wq), the luminosity 
distance for a flat universe is [39] 

D L (z) = cH^il + z) r dz[(l + z) 3 n m + (1 - ft m )(l + *)3(l-H^)]-l/2 ( 19 ) 

Jo 

then the likelihood function depends on Q m and Wq. Fig. [5] shows the likelihood contours 
on (Q m ,wo) plane for GRBs (dark cyan dash-dotted lines), SNe la (blue dashed lines), 
and SNe la + GRBs (red solid lines), respectively. The contours correspond to 1, 2, and 3a 
confidence regions, respectively. The cosmological parameters with the largest likelihood are 
fi m = 0.24+0^4 and w 0 = —0.85lo!5i (l^) with x L™ = 727.32 for 659 degrees of freedom. For 
a prior of Q m = 0.29, we obtain wq = —0.95ld;J|, which is consistent with the cosmological 
constant (i.e., wq = —1) in a la confidence region. 

3.4 w(z) = Wo + W\z/{1 + z) model: Time-dependent Equation of State 

We next examine models in which dark-energy changes with time. As shown above, we 
adopt a simple model in which the dark-energy equation of state can be parameterized by 
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1.0 



Figure 4: The contour confidence levels of (f2 m ,n) in the cardassian expansion model from the 
data for 82 GRBs (z > 1.4) (dark cyan dash-dotted lines), 580 SNe la (blue dashed lines), and 
82 GRBs + 580 SNe la (red solid lines), respectively. The contours correspond to 1, 2, and 3a 
confidence regions. 


[4J Ej 

w(z ) = wo + w\z/ (1 + z) = wo + wi(l — a). (20) 

The ACDM model is recovered when wo = — 1 and w\ = 0. In this dark-energy model, the 
luminosity distance is calculated by 

D L (z) = cH^(l + z) f dz[( 1 + z) 3 n m + (1 - D m )(l + z )3(l+-o+- 1 ) e -3- 1 ,/(l+h]-l/2_ (21) 

Jo 

Fig. [6] shows the constraints on wq versus W\ from la to 3cr confidence regions. The 
dark cyan dash-dotted lines and blue dashed lines represent the constraints from 82 GRBs 
and 580 SNe la, respectively. The red solid contours are obtained from the combination 
of these data. For a prior of Q m = 0.29, we find the best dark-energy parameters set is 
(w 0 ,wi) = (—0.96±o;3g, —0.04l};gg) at the la confidence level with xhm = 727.54/659. This 
result is also consistent with the ACDM model (i.e., Wq = —1 and uy = 0) in the la 
confidence region. 

4 Conclusions and Discussion 

In this paper, we have updated five GRB luminosity relations (ri ag — L, E p — E [so , E p — 
L, E p — £/, trt _ L ) among certain spectral and light-curve features with the latest 139 
GRBs. We find that the five relations indeed exist with the latest GRBs data. To avoid any 
assumption on cosmological models, we obtained the distance moduli of 57 low-z (z < 1.4) 
GRBs by using cubic spline interpolation from the 580 Union2.1 SNe la compiled in Ref. [36]. 
Then, we calibrated the five relations with these 57 low -2 GRBs. In order to constrain 
cosmological models, we extended the five calibrated luminosity relations to high -2 and 
derived the distance moduli of 82 high -2 (z > 1.4) GRBs. 
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Figure 5: Constraints on O m and u>o from la to 3<r confidence regions with dark energy whose 
equation state is constant. The contours are derived from GRBs (dark cyan dash-dotted lines), 
SNe la (blue dashed lines), and SNe la + GRBs (red solid lines), respectively. 


Motivated by these significant updates of the observational data, we considered the joint 
constrains on the Cardassian expansion model and dark energy with 580 Union2.1 SNe la 
sample (z < 1.4) and 82 calibrated GRBs data (1.4 < z < 8.2). In the ACDM cosmology, 
we find that adding 82 high- z GRBs to 580 SNe la significantly improves the constrain on 
Q m — Oa plane. We obtain Q m = 0.27±o;gg and Oa = 0.62+q/q (lcr). For a flat Universe, 
the contours yield = (0.29io!o4) O-^liam)- The transition redshift at which the 

Universe switched from deceleration to acceleration phase is zt = 0.64io!i4 (I 0- )- In the 
Cardassian expansion model, we obtain Q m = 0.24lg;J5 and n = 0.16lo!52 (I 0- )- This result 
is consistent with the ACDM cosmology (n = 0) in the la confidence region. We also fit 
two dark energy models, including the flat constant w model (i.e., w(z) = Wq) and the 
time-dependent w model (i.e., w(z) = wq + w\z/{ 1 + z)). Based on our analysis, it can be 
seen that our Universe at higher redshift up to £ = 8.2 is consistent with the concordance 
model (O m = 0.27, Da = 0.73, w 0 = —l,Wi = 0) within la level. These results suggest that 
time dependence of the dark energy is small even if it exists. 

Since the discoveries of distance indicators of GRBs, these luminosity indicators have 
been used as standard candles for cosmological research at high redshifts. However, the 
dispersion of distance indicators are still large, which restricted the precision of distance 
measurement by GRBs. The large dispersion may be due to that some contamination of 
the GRB sample is unavoidable, and that pure luminosity indicators may never be found 
for these sources. Of course, it could also due to that we simply have not yet identified 
the correct spectral and lightcurve features to use for these luminosity relations. On the 
other hand, it could also due to that we are inevitably suffering from the systematic errors 
and intrinsic scatter associated with the data. In order to estimate distance of GRBs more 
precisely, we should take efforts to investigate possible origins of dispersion of the distance 
indicators, and/or search for more precise distance indicators in the future. 
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Figure 6: Constraints on wq and w\ from lcr to 3a confidence regions with dark energy whose 
equation state is w(z) = wq + w\z/{l + z). The contours are derived from GRBs (dark cyan 
dash-dotted lines), SNe la (blue dashed lines), and SNe la + GRBs (red solid lines), respectively. 
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